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Abstract. We consider the behavior of a modulated wave solution to an § - 
equivariant autonomous system of differential equations under an external forc- 
ing of modulated wave type. The modulation frequency of the forcing is as- 
sumed to bo close to the modulation frequency of the modulated wave solution, 
while the wave frequency of the forcing is supposed to be far from that of the 
modulated wave solution. We describe the domain in the three-dimensional 
control parameter space (of frequencies and amplitude of the forcing) where 
stable locking of the modulation frequencies of the forcing and the modulated 
wave solution occurs. 

Our system is a simplest case scenario for the behavior of self-pulsating 
lasers under the influence of external periodically modulated optical signals. 



1. Introduction. This paper investigates systems of differential equations of the 
type 

^ = fix)+gix)\y\\ (1.1) 
^ = h{x)y + je^"*am, (1.2) 

where a; G M", y G C, the functions /, g : R" ^ M", /i : R" ^ C, and a : M ^ C 
are sufficiently smooth of class C' with some positive integer /. The function a is 
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27r-periodic, and a > 0, /3 > and 7 > are parameters. We assume that for 7 = 
the unperturbed system 

^ = f{x)+g{x)\yf, (1.3) 

f = hix)y, (1.4) 
has an exponentially orbitally stable quasi-periodic solution of modulated wave type 

x{t) = xoiM, y{t) - y^my"'- (1.5) 

Here ao > and /3o > are constants, while : M — ^ R" and yo : R C 
are smooth 27r-periodic functions. We assume that the following nondegeneracy 
condition holds: 

rank diy'^i^) -%o(V') =2. (1.6) 
_ ^y'oW ^yoW 

It is easy to verify that (1.6) is true for all i{j gM. if it is true for one Moreover, 
without loss of generality we assume that ip 1— >■ argyo('0) is periodic, i.e. the curve 
y = yoi'ip) in C does not loop around the origin (otherwise we should replace yoiPot) 
by yo(/3Qi)e*'^^"* and uq by ao — k^^ with an appropriate k G Z). 
It follows from assumption (1.6) that the set 

T2 := {{xom,yo{i^)e"n e x C : ^, ^ e Ti}, 
where Ti = ]R/(27rZ) is the unit circle, is diffeomorphic to a two-dimensional torus. 
Obviously, T2 is invariant with respect to the flow of (1.3)-(1.4), and the solution 
(1.5) lies on 72- 

Roughly speaking, our main result describes the domain in the three-dimensional 
space of the control parameters a, /3 and 7 with |a — aol S> 1 and /3 « /3o such that 
the following holds: For almost any solution {x{t),y{t)) to (1.1)-(1.2), which is at 
a certain moment close to 72, there exists ct G R such that 

\\x{t) - xoiPt + (t)\\ + \\y{t)\ - \yQ{l3t + a)\\ w for large t. 

Let us reformulate our result in a more abstract language as well as in the lan- 
guage of a physical application. 

Abstractly speaking, (1.3)-(1.4) is an autonomous system which is equivariant 
under the Ti-action {x,y) t-^ (x, e*"^?/), ip S Ti, on the phase space. The solution 
(1.5) is a so-called modulated wave solution or relative periodic orbit to the Ti- 
equivariant system (1.3)-(1.4). It is well-known that generically those solutions are 
structurally stable under small perturbations that do not destroy the autonomy 
and the Ti-equivariance of the system. Thus, our results describe the behavior of 
exponentially orbitally stable modulated wave solutions to Ti-equivariant systems 
under external forcings of modulated wave type in the case when the difference 
between the internal and the external modulation frequencies /3 — /3o is small while 
the difference between the internal and the external wave frequencies a — ao is large. 
Note that in [12] related results are described for the case when both differences of 
modulation and wave frequencies are small, and [11] considers the case when the 
internal state as well as the external forcing are not modulated. For an even more 
abstract setting of these results see [4]. 

System (1.1)-(1.2) is a paradigmatic model for the dynamical behavior of self- 
pulsating lasers under the influence of external periodically modulated optical sig- 
nals. For more involved mathematical models see, e.g., [1, 7, 8, 9, 10, 17, 18] and 
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for related experimental results see [6, 15]. In (1.1)~(1.2), the state variables x and 
y describe the electron density and the optical field of the laser, respectively. In 
particular, the absolute value \y\ describes the intensity of the optical field. The 
Ti-equivariance of (1.3)-(1.4) is the result of the invariancc of autonomous opti- 
cal models with respect to shifts of optical phases. The solution (1.5) describes 
a so-called self-pulsating state of the laser in the case when the laser is driven by 
electric currents which are constant in time. In those states the electron density 
and the intensity of the optical field are time periodic with the same frequency. 
Self-pulsating states usually appear as a result of Hopf bifurcations from so-called 
continuous wave states, where the electron density and the intensity of the optical 
field are constant in time. 

The structure of our paper is as follows. The main results are formulated in Sec. 2. 
The proof is splitted into four sections. In Sec. 3 we use averaging transformations [2] 
in order to eliminate the fast oscillating terms with the frequency a. It appears that 
the first non- vanishing terms after the averaging procedure are of order ^'^ jo? . Local 
coordinates in the vicinity of the stable invariant toroidal manifold arc introduced 
in Sec. 4 and then in Sec. 5 the existence of perturbed manifold is proved. The 
global behavior of a system on the perturbed torus is described in Sec. 6. Among 
others, the methods of perturbation theory [13, 14] are used in our analysis. 

2. Main results. In new coordinates x = x,y = re*"^, r, e R, the unperturbed 
system (1.3)-(1.4) has the form 

(iT 

- = f{x)+g{x)r\ (2.1) 
dv 

— = m{x)r, (2.2) 

^=^h{x). (2.3) 

This system has, by assumption, the two-frequency solution 

x{t) =xo{l3f)t), r{t) =ro{(3ot) := \yo{(3ot)\, ip{t) ^ a^t + a.Tgyo{(3ot). 

The subsystem (2.1)-(2.2) does not depend on ip and has an exponentially orbitally 
stable periodic solution x{t) = xoiPot), r{t) = ro{l3ot). The corresponding varia- 
tional system has the following form 

^=A{i^)z, zeM"+i, (2.4) 
dip 

where 

f{xom+g'{xoWyoW 2g{xomroW 

m'{xoi^))ro{^) mMi')) 

We assume that 

the trivial multiplier 1 of the monodromy matrix of linear periodic 
system (2.4) has multiplicity one, and the absolute values of all (2-5) 
other multipliers arc less than 1. 

The adjoint system 



AW ■■= ^ 

Po 
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has a nontrivial periodic solution p{ip) {A^ denotes the transpose of A) , which can 
be normalized such that 



Let us define the function Q : W' 



, , ,s — 1 for all lb- 
rM) J 

X Ti R"+i as follows 

" g{x)\am^ 




(2.6) 



Our first result describes the behavior (under the perturbation by the forcing 
term with 7 > 0) of 72 x M, which is an integral manifold to (1.1)"(1.2) with 7 = 0, 
as well as the dynamics of the system (1.1)"(1.2) on the perturbed manifold. 

Theorem 2.1. Let us assume that the conditions (1.6) and (2.5) are met. 

Then for all /3i < (32 there exist positive constants (jl^,, a*, 5, L and n such that 
for all {a, (3, 7) with 



7 

a > a*, (3i < j3 < 132 and < — < /i* 

a 



(2.7) 



the following holds: 

(i) The system (1.1)-(1.2) has a three-dimensional integral manifold 0Jl(a,/3,7) 
which can he parametrized byil),if,t€M. in the form 

^,2 



X = xoii>) + -^Xi I ip, (f, I3t, at,-,l3, - ] + -^X2 ( ip, ip, l3t, at,— ,{3 



1 



7 



1 



7 



y = ro{ilj)e'^f+'t'^'t'^^ - i-e"''ail3t) 



Yi il},ip,pt,at, -,(3, ■ 



-:^Y2(^,ip,f3t,at,-,(3,l] with 0(V') := r [Qh{xo{0) - c^oH- 
a V a aj fio Jo 



Here Xj 



X U 



and Yj 



X [/ — > C are C smooth, At: -periodic with 



respect to tp and 2TT-periodic with respect to (f, j3t and at and 

U ■.^\{v,|3,^l) (zM? : 0<u < I3i < (3 < 132, Q < ti< ^l 
L a* 

(ii) The dynamics of (1. !)-(!. 2) on 9Jl(a,/3,7) in coordinates ip,(p and t is de- 
termined by a system of the type 

df^ -2 

dt ct^ 



2 4 

\ a/ 



+ ^*2 ( ^p,^,pt,at,-,p,^] I 4>,^,pt,at,-,p,^ \ , (2.8) 



1 



7 



7 



a'^ 



1 



7 



a 



dip 



U, (3t, ^) + ^^2 ^, f3t, at, 

V aJ a'^ \ a a 



7 



1 



+ ^r$3 ip,ip,^t,at, -,(3, - 



where the functions "^1,^1 



X [0,n^) 



and , $j 



X U 



(2.9) 

(J=2,3) 



are C^^'^-smooth, Air-periodic with respect to ip and 2TT-periodic with respect to (p,l3t 
and at. 

(Hi) The integral manifold dJt{a, (3,^) is exponentially attracting (uniformly with 
respect to {a, (3,^) satisfying (2.7)) in the following sense: For any solution (x{t),y{t)) 
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to (1.1)-(1.2) such that dist((a;(to), 2/(^o))i < 5 for certain to G R there is a 
unique solution (V'(i), <p(i)) to (2.8)- (2.9) such that 

2 

x{t) - xoim) - ^Xi{t) - ^X2{t) 

2 

+ yit) ~ z^e^'^'am - ro(V(i))e^^^'*^+'^'^'*"' " W - ^Y2(t) 
< Le-"(*-*«Mist iixito),yito)) ,T2) , t> to, 



< 



whe 



X,(t) -.^ Xj[m,v{t).Pt,at,'^,P,^ , 3 = 1,2, 
Y,{t) := Y, ( ^{t),^{t),pt,at,-,P,^\ , j = 1,2. 



Let us define the function 



and the numbers 



G+ := max G(i6), G_ := min G(V'). 

V'e[0,27r] 'i/'e[0.27r] 



For the sake of simplicity we will suppose that all singular points of G are non- 
degenerate, i.e. 

G"(V') ^ for aU V such that G'(V') = 0. (2.10) 

This implies that the set of singular points of G consists of an even number 27V of 
different points: 

{ij e [0, 27t) : G'(^) = 0} = {^'1, . . . , ^2n}- 
The set of singular values of G will be denoted by 

S:={G{^i),...,G{^2n)}. 

The following two theorems describe the dynamics on D}l{a, f3, 7) in more details. 
In particular, they show that for appropriate parameters (a, /3, 7) there appears 
an even number of two-dimensional integral submanifolds, which determine the 
frequency locking behavior we are interested in. 

Theorem 2.2. Assume that (1.6), (2.5) and (2.10) hold. 

Then for any £ > there exist positive p,* , yU*, and 6 such that for all parameters 
{a,l3,-f) satisfying 

P 



<j<fi^a, (2.11) 



a 

G_ < 4(/3-/3o) <G+, (2.12) 
7 

dist (^^(/3-/3o),5^ >£ (2.13) 



the following statements hold: 
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Figure 1. Graphs of the function G. 



(i) The system (1.1)-(1.2) has an even number of two-dimensional integral man- 
ifolds %{a,l3,j) C an(a,/3,7), j = 1, 27V(a, /3, 7), < N{a,l3,j) < N which 
can be parametrized by tp.t £ W in the form 

X = xo{(it + dj) + -XiALp,l3t,at,-,(i,-] +-X2j['~p,(it,at,-,I3,-^ 
q;\ a a J a \ a a J 

+ ^Yi, (^,l3t,at,-,p,l] +-Y2J (ip,pt,at,-,p,l] , 
a \ a a / a \ a a / 

where dj are constants, and the functions Xkj,Ykj : x — > M are C^^^ -smooth 
and 2iT-periodic with respect to pt and at, and 

V := {{v, /?, pi):G-< - /3o) < G+,^i*i^^ < ^i,,dist{fi\(3 - S) > e} . 

(a) The dynamics of (1.1)-(1.2) on *yij(Q;, /3, 7) in coordinates (p and t is deter- 
mined by an equation of the type 

^^ao + ^^ij (ip,(3t,at,-,p, 1] +^$2j (ip,pt,at,-,(3,l] , 
at a"^ \ a a J a'^ \ a a ) 

where the functions ^kj : M'^ x 1/ — > M are C'""* smooth and 2'K-periodic in (p,pt 
and at. 

(Hi) Any solution {x(t),y(t)) to (1.1)-(1.2) such that dist((a;(to), y(io))j < 5 
for certain € K tends to one of the manifolds Dlj(Q;,/3,7) as t ^ 00. 

Theorem 2.3. Assume that (1.6), (2.5) and (2.10) hold. 

Then for any e > and ei > there exist positive fi* , /i» and S such that for 
all parameters {a,/3,j) satisfying the conditions (2.11)-(2.13) and for any solution 
{x{t),y(t)) of system (1.1)-(1.2) such that dist((a;(io), J/(io))j < S for certain 
to S ]R there exist ct, T e M such that 

\\x{t) - xoiPt + a)\\ + \\y{t)\ - \yoil3t + a)|| < £1 for all t > T. 

The conditions (2.11)-(2.13) from Theorems 2.2 and 2.3 determine the so-called 
locking region, i.e. the set of all triples {a, /?, 7) for which modulation frequency 
locking takes place. These domains are illustrated in the figures 1-4. 
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Figure 3. Cross-sections of the locking region /3 = const. 

In Fig. 1 we show two typical cases of graphs of the function G. In the case (I) 
there exist one positive and one negative local extremum and in the case (II) two 
positive and two negative local extrema, i.e., 

(I) : N^l, S^{G-,G+}, G- <0<G+, 

(II) : iV = 2, S = {G^,Gi,G2,G+}, G_ < d < < G2 < G+. 

In Fig. 2 we show a— const sections of the locking region. In the case (I) this 
section is 

* 2 
(^,7): ^<7<M*a, G^+e<^{l3-l3o)<G+-e 

It is bounded by two straight lines 7 = /i*/a and j ~ fi^a and by two square root 
like curves 

7 = «y ^ ^'^^ G G {G_ + e, G+ - e}. 

In the case (II) the a^const section is bounded by the same two horizontal straight 
lines and by six square root like curves 



-f = a^ ^—J^ with G e {G_ + e, Gi - e, Gi + e, G2 - e, G2 + e, G+ - e}. 
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Finally, in Fig. 3 we show /3=const sections of the locking region in the (1/q!,7) 
plane. Wc consider the parameter a in the region a > a, with sufficiently large 
a* > 



al>^J^. (2.14) 



If (2.14) is satisfied, consider the set of all f3 > /Sq such that 




^ and \^ — — < ^H.. (2-15) 

V G2 — e 

For any fixed a > a*, where a* satisfies (2.14), and for any fixed /3 > /3o with (2.15), 
the line {(a, /3, 7) : 7 e M} crosses the boundary of the locking region in two points 
7 ~ /iia and 7 ~ in case (I) and in four points 7 = /xia, 7 — /i,2a, 7 = fi^a 
and 7 = /i,a in case (II) (see also Fig. 4). Here we denoted 



/3-/3o P-Po P-Po 



3. Averaging. In this section we perform changes of variables with the aim to 
average the nonautonomous terms with fast oscillating arguments at. As the result 
of these transformations, we obtain an equivalent system, where the fast oscillating 
terms have the order of magnitude of jo? and smaller. The principles and details 
of the averaging procedure can be found e.g. in [2]. 
Performing the change of variables 

y = 2/1 - i-e'"^a{pt) 
a 

in (1.1)-(1.2), we obtain the transformed system 

^ = /(xi) +g(a;i)|yip + 45(xi)|a(/3t)p - ^<?(xi)3{2/ie— *a*(/3i)}(3.1) 
at a 

^ - hix,)yi - * Je»* (h{x,)am - /3§(/3t)) , (3.2) 
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where * denotes complex conjugation. In system (3.1)-(3.2), the fast osciUatory 
terms with frequency at are now proportional to j/a. Since the first averaging has 
not produced any nontrivial contributions on the zeroth order, the second averaging 
transformation is necessary: 



7 



which allows eliminating fast oscillating terms of order j/a. 



7^ 7 1 

+ -^ri{x2,y2,at,l3t,-,-), 

^ = /i(.T2)y2 + ^e»*(^2/3/i(x2)^(/3t)-/i2(x2)a(/3t) 

- + ^^{f{^2)+9{x2)\y2?)am) 

- 2^.9(x2)^^5R{2/2e-^"*a*(/3t)} + 4^2(^2, 2/2, e^"*, -, -), 



where the remainder terms ri,r2 are C'~^ smooth functions in all arguments and 
27r-periodic in at and in fit. 

Again, the second transformation has not produced any nontrivial contributions 
of the order l/a. Let us perform the third change of variables 



2^2 



\ dx2 



dxo 



2^g{x2Me'''\h*{x2)ai(3t) + 2/3— m - /i(x2)a(/?i))}. 



dt 



2/2 



2/3 - 



2/i(x-2)/3 — (/3t) - h'{x2)a{l3t) - P'^iPt) 



dh{x2) 

dX2 



{f{x2) + gix2)\y2\^)am 



2^^^.9(x2)5{y2e-"*a*(/3t)}, 
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which transforms the system to the foUowing form: 



dX3_ 

dt 



dt 



= fix3)+gix3)\y3\' 



T 



g{x3)\a{(3t)\' 



7 



7 1 



+ -^'^3 X3,y3,at,l3t, -, - + -^Ti X3,y3,at,l3t, -, - 



7 



7 1 



(3.3) 



= h{x3)yz 

+ [ X3,y3,at,l3t. 



7 1 



re X3,y3,at,f3t 



7 1 



(3.4) 



where the remainder terms r^, ...,re are 27r-periodic in at and in /St, of class C'~^ in 
all variables. The obtained system (3.3)-(3.4) contains a nontrivial contribution of 
the order and all fast oscillatory terms of the orders j/a"^, jc? and smaller. 

The next section proceeds with the analysis of the averaged system (3.3)~(3.4). 



4. Local coordinates. Let us introduce two new parameters 



, c . 

a a 



We assume that /i G {0,fJ.o) and e G (0,£o) with some sufficiently small /io > 0, 
£o > 0. The system (3.3)-(3.4) can be re-written as 



^ = f{x3) + gix3)\y^\' + ^^'g{x3)\am\^ 
dys 



(4.1) 



dt 



+e'^Hr3{x3,y3, 13t, at, /i, e) + e/x^r4(x3, yg, pt, at, fi, e), 
Hx3)y3 + Hr5{x3,y3, pt, at, fi, e) + en'^re{x3,y3, pt, at, fi,e). (4.2) 



After the change of variables 



(4.3) 



in polar coordinates {r,d) the system (4.1)-(4.2) takes the form 
dx 



^ = fix3)+9ix3y +^^^9ix3)\am\^ +£^^'fl+£^^^f2, (4.4) 



^ = m{x3)r + ei^L^f3 + e^fiU, 
dt 

^ = 9/i(x3) + e^i^/s + eV/e, 
at 



(4.5) 



(4.6) 



where fj = fj{x3,r, 9, (3t, at,fi,e),j — 1, . . . , 6, are C'~^-smooth and 27r-pcriodic in 
0,l3t, at functions. Here we assume r > r^, = i min^ l?/o('0)l > 0- 

By substituting z = {x3,r), system (4.4) - (4.6) takes the following form 
dz 

F{z) + ^2^(z, /3t) + e^l^Fl + e^iiF2, 



dt 

d6 
'dt 



h2{z)+efi'^F3 + e^nF4, 



(4.7) 
(4.8) 



where the functions h2, F, and G are defined by h2{z) :^ '^h(x3), 



Fiz) :- 



f{x3) + g{x3)r^ 
m{x3)r 



Giz,pt) 



g{x3)\am\^ 
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Fj = Fj{z, 9, (3t, at, /i, e), j = 1, . . . , 4 are C'~'^-sniooth and 27r-periodic in 9, (3t and 
at functions. The above defined function Q, which is defined for z = (r, x^) € R"+^, 
on the subspace r = 0, i.e. of aU vectors (0,2:3), is just the function Q defined in 
(2.6). Therefore, the use of the same notations should not lead to misunderstanding. 
Equation 

dz , 

has the periodic solution z{t) = ZQ^Sgt) =^ (xo(/3oi)i ''o(/3oO) s-i^d the corresponding 
limit cycle in M"+-'^ is z = zq{iP), ip gTi, i.e., 

—7— = — o , V'eTi. 4.9) 

dip Po 

Let ^i(jp) be the fundamental matrix solution for the variational equation 

dSz _ 1 dFjzo iil;)) 
dtjj dz 

along the periodic solution zq{iP). 

By the Floquet theorem, the fundamental matrix fli{')p) can be represented in 
the form 



-Sz. 



(4.10) 



(4.11) 



where $i(V') is 47r-periodic (n + 1) x (77, + 1) real matrix and Hi is (n + 1) x {n + 1) 
constant real matrix. 

Since dzo{ip)/dip is a periodic solution of (4.10), we can choose 



dzo{ilj) 



$i(^) = 



di' 



<i>(^) 



where ri('0) and $("0) are (n + 1) x n matrices and Hi = diag{0, H} with nxn con- 
stant matrix H. Since the periodic solution zo(V') is orbitally stable, all eigenvalues 
of matrix H have negative real parts. 

Let us find the inverse matrix for $i(?/') : 

$-i(V;) = (ni{iP)e-"^'f'/^<'y^ = e"''f'^f^°n-\ijj). 
Taking into account that 

n'({ip)ni{iij) = I, vgR, (4.12) 

where / is the identity matrix and f2i(V') is the fundamental matrix solution of the 
adjoint system 

we conclude that fl'^^{ip) ~ nj{ip) (see [5]). Accordingly to Floquet theorem 

ni{^l;) = ^i{'il;)e"''f'/P'>. 

It follows from (4.11) and (4.12) that 

nii^p) ^ {ni'wf = ($ri(^))%-^i"^/^o. 

Hence 

i.f(V') = $r'('/'), Hl = -Hi. 
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Since the linear periodic system (4.10) has one nonzero hnearly independent pe- 
riodic solution, the adjoint system (4.13) has also one nonzero linearly independent 
periodic solution. Then 

fii(^) = [p{i^),n{t^)\ , l>i(V') = [p(V'),*(^)' 

where p{ip) is 27r-periodic solution of adjoint system (4.13) and and ^{ip) are 

{n + 1) X n matrices, $(V') is 47r periodic. Taking into account (4.12), we obtain 
that the scalar product in of two vectors dzo{il}) / dip andp(-0) is equal to 1 for 

all V e Ti- 
lt can be verified that 



$i(^). 



(4.14) 



dip dip 
Then (n + 1) x n-matrix ^{ip) satisfies relation 

We introduce new coordinates ip and h instead of z in the neighborhood of the 
periodic solution zq by the formula 

z = zo(V') + $(^)/i, (4.15) 
h\\ < ho with some ho > 0. After substituting (4.15) into (4.7) we 



where h G 
obtain 



fdzoiiP) , d^^)\ dip 



\ dip 



dip 



$(V.) 



dh 



dt ' ' dt 

= F{zo{ip) + <i>iip)h) + fi^g{zo{ip) + Hip)K m 

+eii^Fi{zoiip) + ^{ip)h,e,/3t,at,fi,e)) 

+e^nF2{zo{ip) + ^{ip)h, e, pt, at, e)). 
With regard for (4.9) and (4.14), the relation (4.16) yields 



(4.16) 



fdzoW , d<S>{iP)\ (dip 



V dip 



dip 



/So 



.(,)(f-H. 



= F{zo{iP) + HiP)h) ~ F{zo{iP)) ~ ^^^^$(V)/i 

+^l^g{zo{1p) + <^{ip)h, pt) + e^?Fi{zo{-^) + ^{ip)K 0, f3t, at, 

+eV^2 (2:0(1/0 + ^{'ip)h, 6, pt, at, fi, e). 
Since by our construction det 



(4.17) 



$(V') = det $1 {i') ^ for all iP, the matrix 
h,^iP) 



dip 

dzo{iP) , d<^{iP) 



dip dip 

is invertible for sufficiently small h. Therefore taking into account the expansion 

[A + B)-^ = A-^ - A-^BA-^ + A-^BA-^BA-^ - 
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we obtain for sufBciently small h 



dtp 



*iW 



di) 



h,0 



- p^(^) ■ 








+ 


H2{h,ilj,fi) _ 



where the C'~^-smooth function H{h,ip,fi) — 0{\\h\\) is periodic in ip. 

Hence, the equation (4.17) can be solved with respect to the derivatives dip/dt 
and dh/dt : 



dh 
'dt 

dip 



+ [l>^(^) + H2{h,i:,y.)][F^ + ei?F^ + £^^^2], (4.18) 
: /3o + V(V')e(^o(V'), fit) + ii'H^ih, ^, /i)g(zo(i^), /3i) 

+ [p^{ip) + Hi{h, tp, M)][Ai'Gi + F5 + e/i^Fi + £^^2], (4.19) 



where 



F,{h,i',fi) = F{zo{iP) + m')h) - F{zom - ^^^^<I>(^)/i, 

Gi{h, ip, pt, y) = g(zo(V^) + $(^)/^, /3i) - g{zo{i^), Pt). 
Fj ^Fjizo + <i>ii;)h,e,/3t,at,ti,e), j = l,2. 
We supplement this system with equation (4.8): 
dO 

— = h2{zo{tp) + Hip)h)+efi^F3{zo + ^ip)h,0,l3t,at,fi,e) 
dt 

+e^/iF4(zo + ^iip)h, 6, pit, at, fi, e). 

Using the equality 

2ir 



(4.20) 



^ J h2{zo{ip))di' = ao, 

we replace the angular variable 9 in system (4.18) - (4.20) by (p accordingly to the 
formula 



1 



4, 



[^2(2:0(0) - ao]dC, 



where J"'' is a certain antidcrivative of the function /i2(^o(C)) ~ '^o- 
in the following system 

Hh + fx^Rii + R12 + e/i^i?i3 + e^fiRii, 



As a result we obtain the following system 
dh 



(4.21) 



5 = /3o + ^l^p'^Wg{zoW,^3t) + fi^R2i + R22 + e^L^R23 + £2^i?24(4.22) 
dt 



dip 2 22 

— = aa + fi R31+ R32 + £fJ. R33 + £ 
dt 



(4.23) 
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where functions 

Rii = Rii{h,^P,pt,fx) = [^^i^)+H2{h,4',f^)]5{zoii') + H4')h,m, 
Ri2 = Ri2{hA',l3t,t^) = [^'^ii')+H2{h,^J,fi)]F5^0m\^), 
R21 = R2iih,i',l3t,f^)^Hiih,i;,fi)g{zoW,(3t), 

+ [p^{^)+H,{h,iJ,^^)]G^{h,^P,f3t,^i) = 0{\\h\\) 
R22 = R22{h,i^,^t,^l)^[p'^iiJ) + Hlih,i;,^l)]F5ih,^P,pt,^l)^Oi\\hf), 

i?3i = R3iih,^,l3t,ti)^^[ao~h2izoim{p^WGiMi'),l3t)+R2i), 

Po 

R32 = R32{h,i;,(3t,ii) = h2izo{i;) + <^{i^)h)-h2{zo{i;)) 

-^[h2{zom - ao]R22 = 0{\\h\\) 
Po 

are C'~''-smooth, 47r-periodic in ifj and 27r-periodic in fit. i?i3,i?i4, i?23)-R24, ^33, 
and i?34 are C'~^-smooth functions of {h,tp,if,(3t,at,ij,,e), 47r-periodic in ip and 
27r-periodic in ipl3t, at. 

5. Existence of the perturbed manifold. Using the local coordinates intro- 
duced in the previous section, we investigate here the existence and properties of 
the perturbed manifold. In addition to the circle Ti = R/(27rZ) we will use the 
notation T'^ = R/{AnZ) for the circle of length 47r and Tfe = Ti x • • • x Ti for 

k times 

A:— dimensional torus. 

Lemma 5.1. For € [0,/zo] and e € [0,eo] with sufficiently small /io and Sq, the 
system (4.21)~(4-.23) has an integral manifold 

m^,e = {{h,ij,if,t) : h = u{i},ip,pt,at,^i,e), {ij,ip) eT[xTi,te R}, 

where the function u has the form 

u{ip, ip, pt, at, ^, e) = fi'^uo{tl,', pt, fi) + 

+£/i^ui(V', (fi, f3t, at, fi, e) + e^^U2('0, (f, (it, at, ^, e) (5.1) 

with C^~'^-smooth Air-periodic in -0, 2TT-periodic in ip, fit, at functions uq, ui and U2 
such that ||uj||,^i-4 < Mi, j ~ 0,1,2, where positive constant A/i does not depend 
on a,fi,e. Here \\-\\(ji~4, is the norm of functions from C'^*(T']^ x T3) with fixed 
parameters fj. and e. 

The integral manifold OJl^^e is asymptotically stable in the following sense: there 
exists = i'o(/iO; £0) such that for every initial value {h, ip, ip) at time t with \\h\\ < 
vq, there exists a unique (V'c'ySo) such that 

\\N{t, T, h, Ip, ip) - N{t, T, u{ipo, (po, /3t, ar, fi, e), V'o, 'Po)\\ 

< Le~"(*^"')||(/i, - {u{tpo,Vo,l3T,aT,fi,e),i:o,Vo)\\, t > t, 

where constants L > 1 and k, > do not depend on a,fi,e. N{t,T,h,ij},Lp) is 
the solution of the system (4.21) - (4.23) with an initial value N{T,T,h,^,Lp) = 
{h,il),ip). 



FREQUENCY LOCKING OF MODULATED WAVES 



15 



Proof. Setting (i = /3t, C2 = ctt in the system (4.21) - (4.23), we obtain an au- 
tonomous system 

^ = i?/i + Qi(/i,V',V',Ci,C2,Ai,e), (5.2) 

^ = ^3o + Q2{h,'^p,(p,Cl,C2,^J.,£), (5.3) 

^ ^ ao + Q3{h,ip,(p,CiX2,IJ',£), (5.4) 

where C'~^-smooth functions Qi, Q2 and Q3 are obtained from the right-hand sides 
of (4.21)~(4.23) with regard in fit = Ci,Q:t = C2- The corresponding reduced system 
has the form 

dh dt/j dip dCi dC2 

The eigenvalues of the constant matrix H have negative real parts, hence 

||e^*|l < £6^"°*, t > 0, (5.6) 

where C = const > 1, kq ~ const > 0. 

By introducing new variables C = {ip, f, Cii C2) a-nd new parameters A = (rji, 772, 
773, /i, e) the foUowing system 

^ = ff/i + Qi(/i,C,A), (5.7) 

^=c.o + g(/i,C,A), (5.8) 

coincides with (4.21)-(4.23) if 771 = 772 = ea^^i '?3 = Ci = /^^: C2 = c^^, and 
Qi(/i, C, A) ^ r]iRii + R12 + ?72-Ri.3 + mRu, 
Q = {.Q2,Q3,QA,Qb), 1^0 = (/3o,ao,/3, a), 
Q2 = r]iP^WG{zoW,l3t) + 771^21 + ^22 + '72^23 + r;3i?24, 

= '71-R3I + ^32 + ?72-R33 + ??3^34, Q4 ^ = 0. 

By [13] or [20], for aU parameters A e I\g ~ {A : ||A|| < Aq}, with sufficiently 
small Aq system (5.7)-(5.8) has a unique invariant manifold 

h^woiC^X), CeT'i xTa, Ae/A„, (5.9) 

where u'o(C, A) is bounded Lipschitz in C, A and 'Wo{C^ A) — > uniformly as (771 , 772, 773) - 
0. 

In order to show this, for A G I\„ the mapping T\ : Tp ^ Tp has been used, 
TA^m^ f e-«^Qi(u;(Cr,A),Cr,A)dr, 

where Ct is solution of (5.8) for h = u>(C, A) with initial conditions Co = C- is 
the space of Lipschitz continuous functions w : T[ x R" such that \\w\\c < 

p, Lip w < p, Lip w is Lipschitz constant of w with respect to C- 
Denote 77 = 7/1 + 7/2 + 7/3. Let Mq be a positive constant such that 

< Mo, ||i?^.Ri2|| < \\h\\Hlo, WomisW < Mo, WD^RiiW < Mq, 
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for \\h\\ < po, |A| < Ao, C G T'^ x T3 with some po > 0, Aq > 0. are derivatives of 
order |j| < 1 — 4 with respect to hX,^ (first derivatives of R12 with respect to h 
have estimate ||/i||Afo and higher derivatives have estimate Mq). 

We consider the subset J^r/ao of J^pg which consists of functions w with 5: 
Tjao, Lip(^w < r/ao, where ag is some positive constant. 

For sufficiently small rj, the mapping 

Tx{ri) ■ J^vao ^ ^vao (5.10) 

is well defined. Here A(77) means A = (ryi, r]2, 773, p., e) with 771 + 7^2 + 7?3 = r]. Really, 
for w € -7>;aoi the function Qi has the following estimate 

||Qi(i(;(C, A), C, A)||c < vMo + V^alMo, 
hence, taking into account (5.6), 

£_ 

Let Cr and Cr be two solutions of (5.8) with h = w{(, A), \\w\\c < V^o, Lip(^w < 
Tjao and initial values Cq and Co - Then 

where oi is a positive constant independent on rj. Inequality (5.12) permits to esti- 
mate Lipschitz constant of T{w) : 

rA(„)(«^)(Co')-T^A(,)(w)(Co')ll < 

< [ Ce'''"^ (liphQiLip^w + LipcQi) \\C - CrWdT < 

< — ^ LaoLiphQi + LipqQi) IICo - Coll- (5-13) 



\T^ir,){w)\\c<-{r] + ifal)Mo. (5.11) 



< 



mo 



Ko - ai77 
One can verify that 

LiphQi < V^Ia + SiiaoMo, Lip^Qi < rjMo + ifalMo (5.14) 

if < r;ao. 

There exist positive ao and r^o such that 

imoLiphQi + Lip(;Qi) < 7]ao, ^^^(77 + 77^00) 

Kq — aiT] V / Kq 

for all 1] < riQ. Taking into account (5.14), to this end it suffices 

— — -l>77ao, (1 + 47700 + ao7/) < ao. 

MqL Kq — ai7] 

Hence, mapping (5.10) is well defined for 77 < rjQ. 

Analogously to [20] (Theorem 6.1), we show that the map Tx(ri){w) is a contrac- 
tion of set J-'rjao for all 77 < 771 with some rji < rjQ. The mapping 7a(,,) has unique 
fixed point wq{C,, A) for all A € with 77 < 771. 

Expressions in right-hand sides of (5.11) and (5.13) don't depend on a € [aQ, 00) 
(note, that a is contained explicitly only in equation dC,2/dt = a). Hence, values 
ao and 770 can be chosen independent on a S [ao,oo). By construction, wq satisfies 
||?i'o(C, A)|| < 7700 with positive constant ag independent on a. 

Note that by [20], for sufficiently small A the map 

Ta(w) : C'"^(T'i X T3,]R") ^ x T3,]R") 
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is well defined. 

For proving C'~^ smoothness of integral manifold ifo(C, A) we use the fiber eon- 
traetion theorem [3], p. 127. At first we show that invariant manifold is with 
respect to C. The continuous differentiability with respect to A is proved analogously. 
The smoothness up to C'~^ can be improved inductively. 

Following [3], p. 336, we introduce the set J-^ of all bounded continuous functions 
$ that map T\ x T3 into the set of all n x 4 matrices. Let denote the closed 
ball in with radius p. 



For w e J^7]aoi we consider the map Tl{w, $) : Trjaa x -^7702 ^ -^j? 



jja2 ' 



riK*)(C) = 

dQl{w{Cr,X)Xr,X) 

dh 



<i>(Cr,A) ]WiT,X)dT, 



(5.15) 



where Ct, W{t, A) are solutions of the system 

^=c.o + g(«;(C,A),C,A), 

dW _ dQ{wiC,X)X,X) 9Q(^C,A),C,A) 

dt dc ^ dh ^i^'-^^''''- 

Taking into account the structure of the function Q, we see that 

n5QKC,A),|^^^^ II^OKCAliK^ 



(5.16) 
(5.17) 



dh 



with some positive constant K independent on rj. Choosing rj such that KT]{l+a2) < 
Ko/4 and applying Gronwall's inequality, we obtain 



\\W{t,X)\\ < Afe('^°/4)(*-*°), 

where M is some positive constant. 

Taking into account (5.18) and inequalities 



(5.18) 



< MoT] + Moifal, 



dQi 



dh 



< MoV + 3Mor]ao, 



we get 



\\Tl{w, $)|1 < / £e-«"^(l + rjal + rja^ + 3rjaoa2)vMoMe''°^/^dT < 
ACMMq 



< 



3ko 



-(1 + ?/ag + 7702 + 3ryaoa2)?7. 



There exist 02 > and r]2 such that the last expression is less then 7702 for rj < r]2- 
Hence, the mapping Tl{w, $) is well defined. 
Let us consider the mapping 



(5.19) 



Analogously to [3], p. 337, it can be shown that (5.19) is continuous with respect 
to w. Now we prove that the mapping (5.19) is a fiber contraction. For w G J^r/ao 
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and $i,$2 e J^^ja2 get 

\\Tl{w,^,)-Tl{w,^2)\\ < 



<l 



< /^e-''^(^ll^llllW^2!lll<fl-*2|l + 

+ (11^11 + !l^ll!l$i||)l|W^i-Vl^2||^rfr. (5.20) 
By (5.17), we obtain following estimate for \\Wi — W2II : 

\\W,it,X)-W2it,X)\\< 1^11^11 + ||^||||ci,,||^ \\W,{s,X)-W2is,X)\\ds + 

\\^\\\mm,-'f2\\cds. 

Inserting (5.18) into the second integral and applying the Gronwall's inequality, we 
get 

\\Wiit,X)~W2it,X)\\ < ^^e(''"/2)*||$i-$2||c. (5.21) 

Kq 

Putting (5.21) into (5.20), we obtain 

\\Tl{w,'S>i) - Tl{w,^2)\\ < - <i>2||c, 

where 



<r = 1 + flo H (1 + TjaQ + 7702 + 3770002) 

Kq \ Kq 

We can choose ^ < 1 for sufficiently small 77 hence the mapping (5.19) is a fiber 
contraction. It has unique globally attracting fixed point (wo,u'i). By (5.15), it is 
easy to see that wi(C, A) is bounded uniformly to a e [ao, 00). Repeating [3], p. 296, 
one can show that wq is continuously differentiable and Dwq = wi . 

Taking into account that the invariant manifold (5.9) for 7]i = 772 = 773 = equals 
to zero h = it can be represented as 

h = 77iWo('0,Cl,M,'7l) +%Wi(V;,l^,Cl,C2, A) +773W2(V',¥',Cl,C2, A). 

Note that wq does not depend on C2, '72,7?3, and e, since system (5.7)~(5.8) is 
independent on C,2 for 772 = 773 = e = 0. Taking into account the dependence of 771, 
772, and 773 on /i and e, we obtain that the invariant manifold of (5.2)^(5.5) has the 
following form 

h = ?/(7/;,(^,Ci,C2,/^,e) = M^wo(V;,Ci,Ai) + 

+efi^ui{ip, ip, Ci, C2, e) + eV"2('i/', 'P, (17(2,^, £)■ (5-22) 

Respectively, system (4.21) - (4.23) has integral manifold SUt^^e defined by the func- 
tion m('0, ip, l3t, at, ji, e). 
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Since manifold (5.22) is smooth, it satisfies the following relation 

^(^0 + Q2{u, ip, (fi, Ci, C2, e)) + 'g^(^o + Qsiu, ^, f, Ci, C2, At, e)) 
du du 

+ 1^^ + = Hu + Qi{u,il], 143X1X2, iJ-.e). (5.23) 

CC,i CC,2 

Taking into account this expression and performing the change of variables h = 
h + u('(/', V', Cii C2, ^) in system (5.2)-(5.5), we obtain 

^^{h + Qoik V, V?, Ci, C2, At, e)) /i, (5.24) 

where 

Qo{h,ip,ip, Ci,C2,A^,e)/i = (3i(u + V, Ci,C2,At,e) - ^, <P, Ci, C2, e) 
-^(Q2(w + /i,-0,<^,Ci,C2,At,e) - (52(w, V',<P,Cl,C2,M,£)) 

-^(Q3(w + /i, V', V', Ci, C2, e) - Q3(w, -0, <<5, Ci, C2, Ai, e))- 

The function Qo can be represented as a sum of two terms Qq = Qoi + Q02 such 
that Qoi = C'(ll'^ll) and Q02 = C'(Ai)- Therefore there exist ao > and Ati > such 
that IIQollc < 1^0/ (2/3) for all h and with \\h\\ < oq and fi < fii. Here ap docs not 
depend on fii. Taking into account (5.6) and an estimate of the fundamental solu- 
tion for perturbed linear system [13], we obtain the following estimate for solutions 
of (5.24): 

\\Ht)\\ < ||/i(to)||/:e(-''"+^ll'3»llc)(*-*o) < |j/i(to)||/:e-('"«/2)(t-to)^ (5 25) 

Since u is proportional to fi and ao does not depend on for all small enough fj,i it 
holds hci = ao — supQ<^<^^ ll'ullc > 0. Taking into account that h = h + u, one can 
conclude that for all h with \\h\\ < Iiq the inequality < ao and estimate (5.25) 
hold. 

As result, if a^ < A^i and solution {h{t),ip{t), Lp{t)) of (5.2)-(5.5) satisfies the 
condition ||/7,(to)|| < ha at initial moment of time t ^ to then 

\\h{t)-ui^{t),^{t),pt,at,ti,e)\\ 

< Ce-'^^'-'^^Whito) - u{^{to),^{to), Pto, ato, e)\\ (5.26) 

for all t > to- 

By [13] and [19], the integral manifold S[t^,e is asymptotically stable, i.e. there 
exists ui = vi{fio,£o) such that if p{{h,ip,ip),dyifj,^eiT)) < t^i at time r then there is 
a unique (V'o, 'i^'o) such that 

\\N{t, r, h, ip, ip) - N{t, T, u{ipo, (fio, Pr, ar, /i, e), -00, Vo)l! 

< Le-«(*-^)|| {h, i;, ^) - (m(V'o, ^o, Pt, ar, /i, e), Vo, V'o)!! , t > t,(5.27) 

where constants L > 1 k > don't depend on a,^,e, p{., ■) is the metric in R" x 
T'l X Ti, N{t,T,h,ip,(p) is the solution of the system (4.21)-(4.23) with an initial 
value N{t, r, h, 0, ip) = {h, tp, (p), and 971^, ^(t) is the cross-section of Tl^^e for < = r : 

9Jt^,e(r) = {(M(^,v5,/3T,ar,/i,£),V,V') : e T'^ x TJ. 

Inequalities (5.26) and (5.27) assure the exponential attraction of all solutions 
of (4.21)~(4.23) that start at t = to from a small neighborhood of the unperturbed 
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manifold ft, = to solutions on the perturbed manifold St^^e with the rate of 
attraction, which is independent on fj, £ (0,/^o]i£ G (0,eo];Ct > a,. □ 

Corollary 1. The system (4.7)~(4-8) has the integral manifold 

+e^^l^'^)u2{iJ,v,Pt,at,fi,e),iJ,ip,t) : (V',v)eTi xT;,teM}. 

6. Investigation of the system on the manifold . Substituting the expression 
for the invariant manifold (5.1) into the equations (4.22)-(4.23), we obtain the 
system on the manifold 

^ = /?o + MV(V')e(^o(V'), m + i3t, /i) 

+e^^S'i2(?/', (fi, I3t, at, fi, e)) + e^fj.Si3{i>, V, l^t, at, /i, e), (6.1) 

rl.in 

ao + ^ S21 (ip, (3t, fi) + efi S22 (^, <P, f3t, at, fi, e) + 



dt 

+e'^nS23{ip,'f,^t,at,fi,e), (6.2) 

where C'^'*-smooth functions Sj,j = 1,2,3 are periodic in ip,ip, l3t,at. 
Now we assume that the frequencies /3o and /3 are close to each other 

In the system (6.1)-(6.2), we change the variables according to the formula 

^ = /3i + 

and obtain the following system 

^ = -fi^A + iJ^p'^Wt + ^Pi)g{zo{f3t + ^1), (3t) + /5ii(/3< + 4'i,Pt, fi) 
+e^'^Si2i(3t + ipi,ip, (3t, at, ^, e) + e^^iSiziPt + ■01, <p, Pt, at, ^, e), (6.3) 

^ = QO + yU^5'2l(/3i + 'Ipi,l3t,^) + £/i^S'22(/3t + tpl. If, pt, at, ^i, s) 

dt 

+e'fiS23{f3t + Vi, I3t, at, /i, e). (6.4) 
Performing now the change of variables 

" / r T/ 



^i=i^2 + ^ bMe + V'i)S(^o(e + V'i),0-G(Vi)]dC, 

P Ja 



</5 = <^2 + — / [S2li^ + 1pl,^,tJ-) - S2l{'4>l,fJ')]d^, 

P Jo 



where 

G(^i) ^ r P'^{^ + ^i)g{zo{^ + i^i),OdL 

1 

52i(V'i,Ai) ^ y S'2i(C + 'i/'i,^,M)'^C, 
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the system (6.3)-(6.4) takes the form 



ei?Si2{tp2,V2, fit, at, e) + £^/^5i3(-02, <<52, at, /i, e), (6.5) 



dt 

+e^^S'22(V'2, <<52, at, n, e) + £^/^523(V'2, V'2, at, /i, e), (6.6) 

where the functions in the right hand side are C'^^-smooth and periodic in 9i,fi, (3t, at. 
Together with (6.5)-(6.6) we consider the averaged system 

^^ = -Am^+/.^G(^2), (6.7) 

= ao + M^5*2i(V'2,Ai)- (6.8) 



Denote 



G- := min G(0, G+ := max G(£). 

CG[0,27r] Se[0,27r] 



Then for A = (/3 - ^o)/m^ e [G_, G+] the equation 

A = G(0 

has real solutions. 

Assume that A is a regular value of the map G, i.e. all pre-images ^ = of 
A by G(^) are non-degenerate G'(i?^) 7^ 0. Then the number of pre-images is finite 
and even due to the periodicity of G(^). The signs of every two sequential values 
G'(i?°) and G'(i?^+i) are opposite 

G'i^°k-i)^ak>0, G'(^§,) < 0, k^l,...,N. 

At every interval (i^2fe-i' ''^2fe) t^*^ function G{9) — A is positive and 

min G(e) > A 



for every sufficiently small S. 



Analogously, at every interval ('?2fc' ^^2^+1) ^^'^ function G{9) — A is negative and 



max G(0) < A 

for every sufRciently small S. Due to the periodicity of G{9) we identify ^2N+i with 
and with 1?°^. 

The averaged system (6.7)-(6.8) has 2N one-dimensional invariant manifolds 

no -{(^0,^2) : ^2 eTi}. 

The system on the manifold reduces to 

^ = ao + ^^52i(t9°,M)- 

Manifolds ^2k'^ = li---:^: are exponentially stable and manifolds Il2i^_i,k — 
1, N, are exponentially unstable. 
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Lemma 6.1. There exist /xq > and co > such that for all < fi < fio md 
£ ^ Coy^ the system (6.5)-(6.6) has 2N integral manifolds 

= { (V'2, V2, i) : ^-2 = ^° + Vj{^2,Pt, at, ^, e), ^2 € Ti, t € R} , 

where 



2 



= ^'^VjoiPt, ^i) + evji{(p2,f3t,at,n,e) + —Vj2{(p2,Pt,at,fi,e), 

A' 

with C'^* smooth, periodic in (p2,(3t,at functions Vjk, such that \\vjk\\(ji-4 < M3 
with the constant A/3 independent on a, fi, e. 

The manifolds Il2k, k = I, ■■■,N, are exponentially stable in the following sense: 
there exists Sq such that if \1p20 — '&2k\ — ^0 ^'^'^ ^0 then there exists an unique 

(poi such that for t > to the following inequality holds 

\tp2{t, to, V'20, 'Po) - ^2{t, to, -d^k + V2k{'~Poi,Pto, atQ,fi, e), (poi)\ 

+ \(p2{t,to,1p20,'Po) - (p2{t,to,l32k + V2k{'P01, pto, ato, n, e), ipoi)\ 

< /:2e-^'«^(*-*«' (1^0 - V'oi I + IV'20 - ^^fe - V2k{voi,Pto, ato, ii, e)\) , (6.9) 

where constants C2 > 1 and k,2 > are independent on a, fi, and e. 

The manifolds Il2k-i, k = 1,...,N, are exponentially unstable in the following 
sense: there exists 60 such that if \1p20 — ^2fc-il — '^0 cin-d (po gTi, then there exists 
a unique ipoi such that for t < to the following inequality holds 

|V'2(i, to, -020, <Po) - ■lp2{t, to, l^2fc_l + V2k-l{'P01,Pto, ato, e),ipoi)\ 

+ |(/?2(i,<o, ■020, Vo) - V2{t,to,'02k-i + V2k~i{vai, f3to, ato, fi, e), (poi)\ 

< C^e/-^(*-^o) _ ^oil + |V^2o - ^Ik-i - V2k-l{^ouPto,ato,^l,e)ip,lQ) 

where constants £3 > 1 and K3 > are independent on a, fi, and e. 

Proof. Setting (i — jSt, C2 — at in (6.5)-(6.6) we obtain the following au- 
tonomous system on 4-dimensional torus T4 : 

^ = -A/i2 + ^^G{i,2) + /i'5n(02,Ci, m) 

+£^^S'i2(-02, <P2, Ci, C2, Ai, £) + £^M'5'i3(02, ¥^2, Ci, C2, e), (6.11) 

^ = ao + ^i^S2l{1p2,^J) + /i^S'2i(V'2,Ci,Ai) 

dt 

+efi'^S22{ip2, f2, Ci, C2, A*: e) + eV'5'23(V'2, ¥^2, Ci, C2, (6-12) 

Let us consider a neighborhood of the point ip2 = '&2k where k G {1,...,A^}. 
Neighborhoods of points 02 = ''?2fc-i''^' ~ a-re considered analogously. In 

system (6.11)-(6.13), we change the variables 02 = t^2fc + ^'^'^ introduce the new 
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time T = is^t 



^ = -f3kbi + G2{hi)b\ + /i2 5n(^?°, + 61, Ci, m) 

+e5i2(i9§fc + 61, ^2, ClX2,^^, e) + — ^isl^^afe + 61, ^2, Ci, C2, M, e),(6.14) 



dr /i^ 



+ 521 « + 61, /i) + /^'52l(^°,. + &1, Cl, A^) 



+£^22 (l92fe + 61 , ¥^2 , Cl , C2 , A* , £) + — '§23 (l?2fc + , ^2 , Cl , C2 , e) ,(6. 15) 

A* 

(ir /i^ ' dr /i'^ ' 

where 02(61)5? (G(i?°, + 61) - A) + /3;,5i. 

Extending the system (6.14)-(6.16) by introducing new parameters ryi, 772, 773 
and X we obtain the system 

^ = -/3fc&i + G2(&i)&? + rn ^11(^9!], + 61, Ci, 



+??2.Sl2(^92fc + 61, <P2, Cl, C2, e) + '?3^13(l92fc + ^1 , ¥^2 , Ci , C2 , A* , £)j(6.17) 
X"o + ^21 (t^sfc + 61 , Ai) + ^1 ^21 (i?2fe + 61 , Cl , A*) 



+?72522(i9^fe + fei, (^2, Ci, C2, e) + %^23(^?2fc + 61 , ¥^2 , Ci , C2 , A^, e)j(6.18) 

= X/3, ^ = xa, (6.19) 

which coincides with (6.14)-(6.16) for 771 = /^^, 772 = e, rj^ ~ £^ / fJ-, X = I/a*^- We 
assume that A = (771, 772, 771, /i, e) e Ixo ~ • ll'^ll < -^017 and x ^ Xo with some 
positive Aq and xo- 

Let /3fc G [/3m, /3m] with some constants Pm > /3m > 0. 

We consider the function space 

C'-4(T3 X /a„ X [xo,oo) X [/3m, /3a/]) (6.20) 

of bounded together with their / — 4 derivatives functions w((p2, Ci, C2, X, /3fe) de- 
fined on (v32,Ci,C2) e TTa, A G /ao, X G [xo,oo), /3fc € [/3m,;0M], and mapping 

r(w) = / e^'=^Q4(w(<^2€,Cic,C2C,A,x,^fe),<P2C,Ci5,Cic, A)d^, 



where (54 is the right hand side of (6.17), and (p2^ ~ f2{£., Ci, C2, A), Cif = /3C + Ci, 
C25 = ^^ + (2 is the solution of (6.18)-(6.19) for hi = 7i;(v32, Ci, C2, A, x, /3fc). 

One can verify that the mapping T(w) maps the space (6.20) into itself. 

Analogously to the proof of Lemma 5.1, we apply the fiber contraction theorem 
and show that there exists a unique fixed point 

w = ?7i7;a;i(Ci,x, A) + ?72Wfe2(¥'2, Ci, C2, X, A) + V3Vk3{'P2,CiX2,X, A) (6.21) 

of T{w) in the neighborhood of (0, 0) e C'-'^(T[ x T3) x h„. 

Functions in right-hand side of (6.21) are C'~^ smooth and 27r-periodic in ip2,CiX2: 
such that ll^'fejllc'— * ^ -^2, where positive constant M2 does not depend on A, x 
and /3fc. 
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Respectively, there exist fj,o > and co > such that for all < ^ < /io and 
e < coy^ the system (6.14)-(6.16) possesses the invariant manifold 



bi = /i^Wfci(Ci,A*) +£Wfc2(¥'2,Ci,C2,M,£) H Wfc3(v'2,Ci,C2,/^,e)- (6.22) 

^J■ 

Here we have used the same notations Vki, Vk2, and Vk3 for the functions depending 
on parameters A,/i in (6.21) and the corresponding functions depending on in 
(6.22). 

Therefore the system (6.5)-(6.6) has 2N integral manifolds 

H Vj2{f2,(3t,at,fj.,e): (p2&Ti,teR}. 

The manifolds Tl2ki k = 1, ...,iV, are asymptotically stable [13, 19], i.e. there exists 
~ i^oil^OiCo) such that if p((^/'20j <y'2o)j n2fc(to)) ^ t^o a-t time to then there is a 
unique (^20 such that 

\tp2{t,to,1p20,'P2o) - Tp2{t,tQ,^l^. + V2k{'P20, l3to, ato, fi, e), <f20)\ 
+ \f2{t,tQ,l(;20,'fi20) - ip2{t,to,'02k + V2k{'f 20 , , ato , fl, s) , (f2Q)\ 

< £3e-P^-,3(t-to) (1^20 - <f2o\ + IV'20 - ^2k - «2fc(<^20,/3to,«io,M,£)|)(p.23) 

where t > to, constants £3 > 1 and K3 > are independent on /i,e,a, p{.,.) is 
metric in R x Ti, Ii2k{to) is the cross-section of Ii2k for t = tg. 

Since the function 61 = W2fc(<P2, C17 C2j £) is a smooth invariant manifold of 
(6.14)-(6.16) we obtain 



1^ ( ^ + S2l{&'2k + V2k, + f^'S2l{&'2k + V2k, Cl, M) + 

+£S22{'&°2k +W2fc,V'2,Cl,C2,/i,£) + — S'23(l92A; + ?'2fc , (/'2 , Cl , C2 , 1 
+ ^ 2 + ^7 2 = -Pfc^2fc + G2(V2k)V2k + tl- Sii['d2k+V2kXl,N 

dCi dC2 M 

-|-£5'i2(l?^fe -|-U2fc,</32,Cl,C2,M,£) + — 5'i3(l9^fc +W2fc,</?2,Cl,C2,M,e). (6.24) 



Taking into account (6.24) and making the change of variables 



bl = V2k{(P2,ClX2,^J■,£) +b2 
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in (6.11)-(6.13), we obtain the following system (analogously as in the proof of 
Lemma 5.1) 



dr 

dj£2_ 

dr 



Op 



2 + •§2l(l?2fc + V2k +b2,ll)+ fJ-^S2l{^"2k + V2k + 62, Ci,M) 



+eS22{'d°2k +V2k + &2,<^2,Cl,C2,M,e) H S23i'&2k + ^^2k + &2 , <^2 , Cl , C2 , , 

dCi _ (3 d(2 _ a 
dr fj? ' dr ^fi ' 

with C'^'*-smooth functions Tj of (62, f 2X1X21 A*: £)i periodic in <P2j Ci? C2 and uni- 
formly bounded for 62 from some neighborhood of zero. 

For sufficiently small 62, M^, £, and we can obtain the uniform estimate 



T062 + m'Ti + eT2 + -Ts 



Therefore the following inequality holds 

\h2{t)\ < |fo2(to)|e-^'^(*-*°^ 



(6.25) 



for all 62 (to) such that 162(^0) I ^ &20 with some 620 > 0. Since V2k is a sum of three 
terms proportional to /i^, e, and £^/^ respectively and 620 is independent on these 
parameters, for small enough fi^, e, and e'^/fi, it holds 620 — ||i^2fe||c > > 0. Using 
bi ~ 62 + V2k, one can conclude that for all 61 with |6i| < 60 the inequality I62I < 620 
and estimate (6.25) holds. 

As a result, if < /i < ^0 and s < c^^fji and solution (7/;2(t), 'y32(t)) of the system 
(6.5) - (6.6) satisfies the condition |'!/'2(to) ~ '^\k\ — ^0 at initial moment of time to 
then 

\tp2{t) - i?2fe - V2k{v2{t),l3t, at,n,e)\ 

< \Mto) - ^2k - v2k{Mto),Pto,ato, f^,e)\e-^'^^'-*'>'> (6.26) 

for all t > to. 

Inequalities (6.23) and (6.26) assure the exponential attraction of all solutions 
of (6.5)^(6.6) that start at t = to from a small neighborhood of the unperturbed 
manifold 1^2 = '&2k solutions of the perturbed manifold Il2k according to the 
estimation (6.23). 

Considering the system (6.5)-(6.6) in the neighborhood of the manifolds n2A:_i, 
k ^ I, . . . , N, we obtain similarly that these manifolds are exponentially unstable 
according to (6.10). 

Corollary 2. The system (6.1)-(6.2) has 2N integral manifolds 

V^ = {{f3t + ^° + ij{^,f3t,at,fi,e),'f,t): ^eTi,teM}, (6.27) 
where the C''^'^ -smooth function 



Vj = fi Vjoil3t,fi) + evji{tp,pt,at,n,e) H Vj2{ip, I3t,at, ii,e), 
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is periodic in (p, (3t, and at. On the manifolds ( 6.27), the system (6.1)-(6.2) reduces 
to 



^ = ao + fi^S2iWt + + Vj,Pt, ^) + e^?S22{Pt + 7?° + ij,ip, (3t, at, fi, e) 

+ e^fiS23il3t + d° + Vj,ip,l3t,at,n,e). (6.28) 

The manifolds corresponding to j = 2k, k — I, N, are exponentially stable for t — > 
+0O and the manifolds corresponding to j ~ 2k — l,k = 1,...,N, are exponentially 
stable for i — > — oo. 

Lemma 6.2. There exist > and cq > such that for all < fi < iJ,o and 
£ ^ coyT^ the solutions of (6.5)-(6.6) have the following properties: 

(i) if a solution {il'2{t),tp2{t)) at a certain time t ^ to has the value ?/'2(io) = 
'^2A;-i + '^1 then it reaches the value '02(^0 + T) = 'd2k ^ ^ after a finite time interval 
of the length T = T{d, n, e); 

(a) if a solution {'tjj2{t), ip2{t)) at a certain time t = to has the value 'ip2{to) = 
'&2k+i ^ '^i then it reaches the value ^'2(^0 + T) = 'd2k + ^ after a finite time interval 
of the length T ~ T{S,^,e). (Here we identify '&2N+1 with i?]* and -do with 1^2^). 

Proof. Let us consider the interval (^?2fe-ii '^^/c)- The intervals (i?2fei^2fe+i) 
be considered similarly. Denote 

m = min (G(£) - A) > 0. 

The right-hand side of (6.5) can be estimated as follows 
dM-t) 



dt 



M^(G(V2) - A) + + sfi^Si2 + s^f^Su >fi'^{m- mo) , 



^2 ~ 



where mo = mo{fJ.o,co) = supyfj,^Sii+eSi2 + ^Si3j . By choosing sufficiently 
small /^o and co , one can obtain toq < ™- Hence 

'02(i) > -02(^0) + M^("i - mo){t - to), 

, . ^ Mt)-Mto) ^ ^2k-^°2k-i ^ 27r 

t-to< < —, r < —, 7 = T[5, fi, e). 

yU^(m — mo) fj, (m ^ mo) /.i^(?7i — moj 

Proof of Theorem 2.1. Theorem 2.1 follows from Lemma 5.f and the following 
chain of coordinate changes: averaging transformations from section 3, (4.3), and 
the local coordinates (4.15) in the neighborhood of the invariant manifold T2. 

Proof of Theorem 2.2. In Lemma 6.1, the existence and local stability prop- 
erties of the integral manifolds IIj, j = 1,...,2N have been proved. The integral 
manifolds correspond to the manifolds Hj after the averaging and transforma- 
tions (4.3) and (4.15). 

It has been proved in Theorem 2.1 that all solutions from some neighborhood of 
the torus 72 are approaching the perturbed integral manifold D)l{a, (3, 7). Therefore, 
for the proof of the statement 2 of Theorem 2.2 it is enough to show that the 
solutions on this manifold are approaching the solutions on one of the manifolds 

Let us fix any positive ei. For the set S of singular values of G we define two 
following sets: 

B{ei) = {g e [G_,G+]; dist{g,S) > ej, 
A{ei) ^{9e [0, 27r] : G{e) e B{ei)}. 
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Taking into account that the sets B{ei) and A{ei) are compact one can prove 
that there exists a positive constant such that 

>^ for all eeA{ei). (6.29) 

clu 

Let us consider the system (6.5)-(6.6), which describes the dynamics on the 
manifold /3,7). For any a, 7 and (3 satisfying (2.12) and (2.13) there exists a 

finite number of points j = 1, 2N, (solutions of the equation /?— /3o = IJ.^G{9)), 
which define the integral manifolds n°, j — 1, 2N, of the averaged system (6.7) - 
(6.8). Note that number N depends on the parameters a, 7 and /3. 

By Lemma 6.1, for fixed A G S(£i), there exist > and cq > such that 
for all < /i < /^o and e < coy/Jl the system (6.5)-(6.6) has 2N integral manifolds 
Ilj. Due to the uniform estimate (6.29), it follows from the proof of Lemma 6.1 
that constants ^0 > and cq > can be chosen the same for all A e B{ei), and 
therefore for all a, 7 and j3 satisfying (2.12) and (2.13). 

All the manifolds Il2k are asymptotically stable in the sense of the formula (6.9) 
and the manifolds n2fe_i, 1 < fc < iV are asymptotically unstable in the sense of 
the formula (6.10). Therefore, if |V'20 — ^^2^-1 1 ^ ^0 and (^^20, ¥'20) ^ ^2k-i then 

Mt)] > A'2e^'''^(*-*°)|62(io)|, t > to, (6.30) 

where K2 > 1, K2 > are some constants and &2(i) = '^2it,to,ip20,V20i£, fJ-) — 
ip2{t, to, -ff^k-i + ev2k-iiv'2o,f3to, ato, e, n),^2o, e, m)- 

It follows from (6.30) that on a finite time interval T depending on values V'20 
and yu, e the solution (V'2(i), f2{t)) of (6.5)-(6.6), whose initial value ("02(^0), ^2(^0)) 
for t — to does not belong to the manifold Il2k-i, i.e. 

if^2{to) 7^ '&2k-i + v2k-i{ip2{to), Pto,ato,e, ^j.), 

and \tp2{to) — ^2k-i\ < '^Oi reaches the boundary of (5o-neighborhood of i^^fc-i' nrore 
exactly, values ip2{,ti) = -ff^k-i ^ '^0 or ^'2(^2) = '&2k~i + '^o- 

Then, by Lemma 6.2, on a finite time interval, this solution reaches (Jg-neighborhood 
of point i?2fc O'^j respectively, ^o-neighborhood of point i?2fc+2' '^here So is defined 
from Lemma 6.1. 

Next, by Lemma 6.1, as t further increases, the solution is attracted to one of 
the stable integral manifolds Il2k or 112^+2. 

As a result, solutions (ipit), (p{t)) of the system (6.1) - (6.2) that, at initial point 
t = to do not belong to the unstable integral manifolds V2k-i, k = 1, ...,iV, i.e., 

V'(^o) 7^'d2k~i + + Tj2k-i{'p{to), Pto,Oito, n,e), 

are attracted for t > to to solutions {ip{t), <f>{t)) on one of the stable integral mani- 
folds V2k 

^{t) = (3t + i9°2k + V2k{'p{t), I3t, at, fi, e), 

(p{t) is a solution of system (6.28) for j = 2k, 

so that 

m)-m\+\^{t)-m\ < C2e-^^'-^-^'-^^ (mr) - ^pin + \^{t) - ^{t)\) ,t>T, 

for some T = T{ip{to), e) and some £2 > 1- 
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If a solution {ip{t), ip{t)) of (6.1)-(6.2) at the initial point t ~ to belongs to one 
of integral manifolds V2k+i then this solution has the following form 

^{t) = /3t + + V2k+i{v{t).pt, at, fi, e), 

(p{t) is a solution of system (6.28) for j = 2k + 1. 

Using the last formulas and Lemma 5.1, we conclude that any solution {h{t), ip{t), ip{t)) 
of (4.21) - (4.23) that starts from the I'o-neighborhood of the integral manifold Ti 
is attracted to one of the solutions (h{t),ip{t), (p{t)) on the integral manifold SOt^^e 
such that 

h{t) = u{tp{t),(f{t),l3t, at, fi, e), 

■^{t) = l3t + ^° + ij{^{t),l3t,at,fi,e), 

(p{t) is a solution of system (6.28) 

with some j, 1 < j < 2N. More exactly, there exist constants L > 1,k > and 
T = r(/i(to), ipito), ^{to)) > to such that for i > T : 

\h{t) - uiiit), ^{t),f3t, at, fi,e)\ + m) - ^{t)\ + \^{t) - ^{t)\ < 
< Le-«(*-^) {\hiT) - u{i^{T),^{T), (iT, aT, ^l,e)\ + 

+|^(T)-^(r)| + |^(r)-^(r)f 



Proof of Theorem 2.3. Under the conditions of Theorem 2.3, the conditions of 
Theorem 2.2 are satisfied. Therefore, every solution {x{t),y{t)) of the system (1.1)- 
(1.2) that at a certain moment of time to belongs to a 5-neighborhood of the torus 
72 tends to some solution on one of the integral manifolds 91j (a, /3, 7), j ~ 1, 2N. 
Hence, for any e > the following inequality holds 



x(i)-xo(/3i + 



|y(t)|-|yo(/3i + <)| 



< e 



with some 1 < j < N for all moments of time starting from T(x{to),y(to)). 

Acknowledgments. The authors gratefully acknowledge the scientific cooperation 
and the helpful discussions with K. R. Schneider over many years which lead to the 
two preliminary versions [16, 14] of the present research. 

REFERENCES 

[1] U. Bandelow, L. Recke and B. Sandstede, Frequency regions for forced locking of self-pulsating 

multi-section DFB lasers, Opt. Commun., 147 (1998), 212-218. 
[2] N. N. Bogoliubov and Yu. A. Mitropolskii, "Asymptotic Method in the Theory of Nonlinear 

Oscillations," Gordon and Breach, New York, 1961. 
[3] C. Chicone, "Ordinary Differential Equations with Applications," 2"^* edition, Springer- 

Vcrlag, New York, 2006. 
[4] D. Chillingworth, Generic multiparameter bifurcation from a manifold, Dyn. Stab. Syst., 15 

(2000), 101-137. 

[5] B. P. Demidovich, "Lectures on Stability Theory," Nauka, Moscow, 1967. 

[6] U. Feiste, D. J. As and A. Erhardt, 18 GHz all-optical frequency locking and clock recovery 

using a self-pulsating two-section laser, IEEE Photon. Technol. Lett., 6 (1994), 106-108. 
[7] M. Lichtner, M. Radziunas and L. Recke, Well-posedness, smooth dependence and center 

manifold reduction for a semilinear hyperbolic system from laser dynamics. Math. Methods 

Appl. Sci., 30 (2007), 931-960. 
[8] M. Nizette, T. Erneux, A. Gavrielides and V. Kovanis, Stability and bifurcations of periodically 

modulated, optically injected laser diodes, Phys. Rev. E, 63 (2001), Paper number 026212. 



FREQUENCY LOCKING OF MODULATED WAVES 



29 



[9] D. Peterhof and B. Sandstede, All-optical clock recovery using multisection distributed- 

feedback lasers, J. Nonlinear Sci., 9 (1999), 575-613. 
[10] M. Radziunas, Numerical bifurcation analysis of the traveling wave model of multisection 

semiconductor lasers, Physica D, 213 (2006), 98-112. 
[11] L. Rccke, Forced frequency locking of rotating waves, Ukrain. Math. J, 50 (1998), 94-101. 
[12] L. Rccke and D. Peterhof, Abstract forced symmetry breaking and forced frequency locking of 

modulated waves, J. Differential Equations 144 (1998), 233-262. 
[13] A. M. Samoilenko, "Elements of the Mathematical Theory of Multi-Frequency Oscillations," 

Kluwer Acad. Publ., 1991. 
[14] A. M. Samoilenko and L. Recke, Conditions for synchronization of one oscillation system, 

Ukrain. Math. J. 57 (2005), 1089-1119. 
[15] B. Sartorius, C. Bornholdt, O. Brox, H.J. Ehrko, D. Hoffmann, R. Ludwig and M. Mohrle, 

All-optical clock recovery module based on self-pulsating DFB laser. Electronics Letters 34 

(1998), 1664-1665. 

[16] K. R. Schneider, Entrainment of modulation frequency: a case study. Int. J. Bifurc. Chaos 

Appl. Sci. Eng. 15 (2005), 3579-3588. 
[17] J. Sieber, Numerical bifurcation analysis for multisection semiconductor lasers, SIAM J. 

Appl. Dyn. Syst. 1 (2002), 248-270. 
[18] S. Wieczorek, B. Krauskopf, T. B. Simpson and D. Lenstra, The dynamical complexity of 

optically injected semiconductor lasers, Phys. Rep. 416 (2005), 1-128. 
[19] Y. F. Yi, Stability of integral manifold and orbital attraction of quasi-periodic motion, J. 

Differential Equation 103 (1993), 278-322. 
[20] Y. F. Yi, A generalized integral manifold theorem, J. Differential Equation 102 (1993), 153- 

187. 

Received xxxx 20xx; revised xxxx 20xx. 

E-mail address: reckeOmath.hu-berlin.de 
E-mail address: samaimath. kiev.ua 
E-mail address: teplinskyaimath.kiev.ua 
E-mail address: vitkSiinath. kiev.ua 
E-mail address: yancluikamath.hu-berlin.de 



